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^ i Abstract 

O 

We investigate the escape rate of the biaxial nanospin particle with a magnetic 



ctf 



O 
O 



field applied along the easy axis. The model studied here is described by 
the Hamiltonian H = -AS 2 Z - BS 2 - HS Z , (A > B > 0). By reducing this 



Hamiltonian to a particle one, we derive, for the first time, an effective particle 



potential for this model and find an analytical form of the phase boundary 
line between first- and second-order transitions, from which a complete phase 

X" 

diagram can be obtained. We also derive an analytical form of the crossover 

temperature as a function of the applied field at the phase boundary. 
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Recently the quantum-classical phase transition of the escape rate in nanospin sys- 
tem has been studied intensively. One of the main issues in this subject is to determine 
whether the transition is first-order or second-order. In this regards, for the uniaxial spin 
system such as high-spin molecular magnet Mn^Ac ||, two models have been investigated 
comprehensively : one with a transverse field || and the other with an arbitrarily directed 
field 1, described by the Hamiltonians H = -DS 2 - H X S X and H = -DS 2 - H X S X - H Z S Z 
respectively. In the first case, by using the method of particle mapping and the Landau 
theory of phase transition, Chudnovsky and Garanin have shown that the transition order 
changes from first to second when the field parameter h x = H X /(2SD) is 0.25. In the case 
of model with arbitrarily directed field Garanin et al. have obtained the phase boundary 
line h xc = h x (h z ) to show the whole phase diagram in which h xc (0) = 0.25 in the unbiased 
case and h xc ~ (1 — h z ) 3 ^ 2 in the strongly-biased limit. 

For biaxial spin system such as iron cluster Fes @ Liang et al. considered a model 
without an applied field, H = K(S 2 + XS 2 ), (0 < A < 1) ||. Using the coherent spin state 
representation they have shown that the coordinate dependent effective mass leads to the 
first-order transition and the change between the first- and second-order transitions occurs 
at the value A = 0.5 . The biaxial spin model with transverse field, Ti = K(S 2 + \S 2 ) —H y S y , 
has also been investigated by Lee et al. who demonstrated that various types of combinations 
of the first- and second-order transitions are possible depending on A and H y /K\S f7|. 

In this paper we study the phase transition of the escape rate of the biaxial spin system 
with a longitudinal field. We will first derive an effective particle potential by mapping the 
spin Hamiltonian onto particle one, which is the first derivation for the present model. Then, 
with the help of the recently developed method for the criterion of the transition order |||| , 
we find an analytical form of the phase boundary line from which a complete phase diagram 
for the order of the phase transition is obtained. 

Consider a nanospin particle with an applied field H along the easy axis. If the spin 
particle is a biaxial spin system with XOZ easy plane anisotropy and the easy Z-axis in the 
XZ-plane the Hamiltonian can be described by 



H = -AS 2 Z - BS 2 X - HS Z (1) 

where the anisotropy constants satisfy A > B > 0. Our model is equivalent to H. — 
K(S 2 Z + XS 2 ) - HS X ,(X < 1) if we set A = K, B = (1 - \)K. In the following, for 
convenience, we introduce dimensionless anisotropy parameter b = B/A and field parameter 
a = H/SA, (0 < a < 2) where S is the spin number. For iron cluster Fes m Ref.5, 
S = 10, A — 0.316 K, and B = 0.092 K. We can reduce this spin problem to a particle 
moving in a potential ]TD|. The equivalent Schrodinger-like equation is derived as 
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is the particle wave function, and V(x) is the effective particle potential given by 
V(x) a 2 S 2 cn 2 x - 2abS(2S + l)snx - 465(5 + 1) 



(4) 
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in which snx, cnx, and dnx are the Jacobian Elliptic functions with modulus k 2 = 1 — b. 
This potential is shown in Fig.l. The local minimum represents a metastable state of the 
spin system described by the Hamiltonian (1). The potential preserves the symmetry of a 
rotation from — Z to +Z-axis. Thus, the escape from the local minimum — x m to the global 
minimum x m corresponds to the inversion of the spin magnetization vector. For large spin 
such as S(S + 1) ~ S 2 = (5 + 1/2) 2 it has a maximum at point x = sn _1 [— a/2(l — b)]. For 
a given value of b the barrier height decreases as a increases and vanishes at a = 2(1 — b) 
which corresponds to the coercive field. 

The type of the phase transition is determined by the behavior of the Euclidean time 
oscillation period r(E), where E is the energy, near the bottom of the Euclidean potential 



which corresponds to the top of the potential barrier |]]. When t(E) decreases monotonically 
as E approaches the barrier top the transition from quantum to classical regime is second- 
order. If, however, r(E) is not a monotonous function of energy a first-order quantum- 
classical transition takes place. One can also argue that the condition for the first-order 
transition can be obtained by looking at the behavior of the oscillation period in Euclidean 
time as a function of oscillation amplitude near the barrier top. In this case the first-order 
transition appears when the amplitude dependent period r(a), where a is the amplitude, 
is smaller than the zero amplitude period r(0) which corresponds to the solution near the 
position of the sphaleron solution ||. The important feature of these approaches relies on the 
shape of the potential near the top of the barrier. Thus, parameterizing the amplitude as the 
coefficient of the perturbation expansion a sufficient condition for the first-order transition 
can be derived 0||. Below, we show that the latter approach leads to an analytical form of 
phase boundary between first- and second-order transitions for the present model. 
Expanding the potential V(x) near x = x up to fourth-order we obtain 

Viz + x Q ) = a 1 (x )z 2 + a 2 (x )z 3 + a 3 (x )z' i } (5) 

where z = x — xo,ai(xo) = V"(xo)/2(< 0), 02(^0) = V'"(xo)/6(> 0), and a 3 (xo) = 
V""(xq)/24. Following references 8 and 9 the criterion of the first-order transition can 
be obtained from the condition r(a) — r(0) < 0, which is expressed by 

_15a!^ Q (6) 
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By equating both sides we can find an equation of the phase boundary line. For large spin 

number this is obtained to be 



« C = 2(1-6 C )^/^, (7) 

where a c and b c are the critical values of the field and anisotropy parameters on the phase 
boundary, respectively. In Fig. 2 we have plotted b c as a function of S c = 1 — a c /2. This 
picture displays a complete phase diagram in the whole range of the applied field. From the 



condition in Eq.(6) the first-order exists below the line. From the picture we immediately 
see b c = 0.5 in the unbiased case S c — 1. This is the same as the previously obtained result 
H and confirms that the present approach is correct. In the strongly biased case, 5 C — > 
and b c — > since a c — > 2. Thus, in this limit, by expanding the Eq.(7) for small b c we find 
the linear behavior b c ~ 0AS c . 

For the critical temperature at the boundary between the first- and second-order transi- 



tions we use the formula 2ttT c = J—V"(xo)/m. For the present potential this becomes 

T c 2y/Sb c ll-b c 



(8) 
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In Fig. 3 we have plotted T C /SA vs. 5 C graph, where the relation in Eq.(7) has been used. 
In the unbiased case Eq.(8) gives T C /SA = 0.318 which coincides with the value calculated 
from the Eq.(13) in Ref.6. In the strongly biased limit it represents linear behavior T C /SA ~ 
0.442<5 C . 

Very recently the same spin system has also been considered ||11| , but with a slightly 



different model : H = —DS^ + BS% — HS Z , (D, B > 0). Comparing this with our model we 
realize D = XA. By using a perturbative approach with respect to b = B/D they obtained 
a phase diagram in the whole range of field from which the linear dependences of b c and 
T C /SD on S c in the strongly biased case can be found. Since the perturbation parameter is 
b the validity of this approach is limited in the range of small values of b which requires the 
strong bias limit to change the order of transition. On the other hand, there is no restriction 
to b in our approach. We believe that the present approach is more rigorous and the results 
are improved. 

Finally, we comment on the experimental observation of the results. Using the anisotropy 
constants given in Ref.5 we have b c = 0.29, and thus a c = 0.81 from Eq.(7) for which the 
critical field is estimated to be 1.9 T, and the coercive field is 3.4 T. From Eq.(8) the 
transition temperature on the phase boundary can also be calculated, and we find T c = 0.79 
K. 

To summarize we have investigated the phase transition of the escape rate from 



metastable state in nanospin system with a magnetic field applied along the easy axis. 
By using the particle mapping we derived an effective particle potential. We obtained an 
analytical form for the equation of the phase boundary line between the first- and second- 
order transitions and thus a complete phase diagram. In the strongly biased case we found 
a linear dependence of b c , the dimensionless anisotropy parameter on the applied field. We 
also obtained a diagram for the crossover temperature as a function of the applied field. 
It also shows a linear relation. The results obtained here can be used as a guide for the 
experimental observation. 
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FIGURES 

FIG. 1. Effective particle potential 

FIG. 2. The phase boundary line b c vs. 5 C . Solid line: anaytical result for large spin number, 
dashed line: straight line in the strongly biased case. 

FIG. 3. Critical temperature T C /SA as a function of S c . Respective lines represent the same as 
Figure 2. 
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